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ABSTRACT 

Modelling the build-up of halos is important for linking the formation of galaxies with 
cosmological models. A simple model of halo growth is provided by Press-Schechter 
(PS) theory, where the initial field of density fluctuations is smoothed using spherically 
symmetric filters centred on a given position to obtain information about the likelihood 
of later collapse on varying scales. In this paper the predicted halo mass growth is 
compared for three filter shapes: Gaussian, top-hat and sharp fc-space. Preliminary 
work is also presented analysing the build-up of halos within numerical simulations 
using a fricnds-of-friends group finder. The best-fit to the simulation mass function 
was obtained using PS theory with a top-hat filter. By comparing both the backwards 
conditional mass function, which gives the distribution of halo progenitors, and the 
distribution of halo mergers in time, the build-up of halos in the simulations is shown 
to be better fitted by PS theory with a sharp fc-space filter. This strengthens previous 
work, which also found the build-up of halos in simulations to be well matched to 
PS theory with a sharp fc-space filter by providing a direct comparison of different 
filters and by extending the statistical tools used to analyse halo mass growth. The 
usefulness of this work is illustrated by showing that the cosmological evolution in the 
proportion of halos that have undergone recent merger is predicted to be independent 
of mass and power spectrum and to only depend upon cosmology. Recent results from 
observations of field galaxies are shown to match the evolution expected, but are not 
sufficiently accurate to usefully distinguish between cosmological parameters. 

Key words: galaxies:formation, galaxies: halos, galaxies:evolution, cosmology: theory, 
cosmology: dark matter 
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1 INTRODUCTION 

Deep surveys are beginning to probe the morphological evo- 
lution of non-active field galaxies out to redshifts where 
cosmological changes in the population as a whole become 
apparent. The evolution of such galaxies is expected to be 
driven by the hierarchical build-up of the dark matter halos 
in which the galaxies reside, with mergers between halos of 
approximately equal size being particularly important. Such 



mergers can lead to tne disruption ot galaxy structure pro- 
ducing observable signatures. C6nsequently, it is important 
to be able to model the growth of halos in sufficient detail 
to determine the relative frequency of mergers. 



Press-Schechter (PS) theory (Press & Schechter 1974) 



provides a simple analytic model for the build-up of ha- 
los. In this theory, the initial field of density fluctuations 
is smoothed on varying scales around a given position to 
obtain information about the probability of later collapse. 
The distribution of overdensities as a function of filter size 
is called a trajectory. If filter size and halo mass are linked, 
each trajectory contains information about the mass of halo 



within which a test-particle associated with that trajectory 
and position in the linear density field resides. Linking per- 
turbation amplitude with the epoch at which fluctuations 
undergo non-linear collapse, and assuming that collapse oc- 
curs on the largest scale for which there is sufficient excess 
overdensity provides a complete model for halo mass growth. 
The standard use of PS theory is to calculate the mass 
function f(M\t) dM, t he distribution of mass in i solated ha- 
los at a given epoch (Press &, Schechter 1974; Peacock & 



Heavens 1990; Bond et al. 1991). For the sharp fc-space fil- 
ter, it is possible to obtain an analytic formula, although 
this is not true for more general filters. A numerical method 
for generating trajectories for ge neral spherically symmetric 
filters was given by Bond et al. (1991) and is summarized in 
Section @. The effect of the choice of filter on the resulting 
mass function is considered in Section 0, where the top-hat 
filter is shown to provide the best match to recent N-body 
simulation results. The distribution of times at which halos 
of a giv en mass are created c an also be determined by PS 
theory (Percival fc Miller 1999), and is discussed for different 



filters in Section 
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It is possible to use PS theory to follow the build-up 
of halos that have mass M at time t by considering the be- 
haviour of trajectories that pass through this location. This 
is often called 'extended PS theory' and leads to cond itional 
anal ogues of the mass function and halo creation rate ( Bower 
199l|; glond et al. 199ll Lacey & Cole 19931 Percival & Miller 
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199E). Extended PS theory is often used to estimate halo 
merger trees, which are i n turn used in semi-analyti c mod- 
els of galaxy fo r mation (|Ka.nffmarm fc White 1993 Shetb] 
& Lemi on 1999J ; fcomerville" et al. 200d| ). The distribution of 
progenitor masses at an earlier epoch is termed the back- 
wards conditional mass function, while the distribution of 
halo masses in which a particle might find itself at a later 
time is te rmed the forward conditional mass function. Lacey 
& Cole (1994) showed that the forward conditional mass 
function from simulations with power law spectra appears 
to be well fitted by a variant of PS theory. In this paper this 
work is extended to determine the effect of filter choice on 
the backwards conditional mass function and the predicted 
frequency of merger events as a mechanism for the build-up 
of dark matter halos. The backwards mass function is anal- 
ysed in Section ]A Starting with the premise that mergers 
between halos correspond to jumps in the mass associated 
with a single trajectory in PS theory, the relative frequency 
of mergers are compared for different filters in Section pi. 

The usefulness of this analysis is demonstrated in Sec- 
tion H, where the proportion of halos that have undergone 
recent merger is calculated. This is linked to the field galaxy 
merger rate, and recent observational results are shown to 
be consistent with a variety of cosmological models. Further 
observations could, in principle, lead to estimates of Q m and 

i Li] . 



2 THE PRINCIPLES OF PRESS-SCHECHTER 
THEORY 

Press-Schechter (PS) theory is based on two simple ideas. 
First, the field of initial density fluctuations is filtered to 
determine the amplitude of the perturbations on different 
scales. Second, the perturbation amplitudes are linked to 
their collapse time or lack of it, thereby providing a com- 
plete statistical model for the build-up of structure in the 
Universe. Three filters are compared in this paper, defined 
in k- space by 

Sharp fc-space 

W ks (k;R ks )= $(l-kR ks ) (1) 

Top-Hat 

W th (k;R th )= 3 s Mfc^-(^)cQ S (fc^) (2) 



(kRth) 



Gaussian 
W 9 (k;R g ) 



exp 



-R 2 g k 2 



(3) 



where $(x) is the Heaviside step function, $(x) = 1 for x > 
and -d(x) — otherwise. The mass variance, a M inside a 
volume V whose shape and size are dependent on the filter 
is given by 



where P(k) is the power spectrum of the initial density per- 
turbations, and is often simply referred to as the power spec- 
trum. Unless stated other wise, the standard CDM power 
spectrunr of Bond et al. (1991) has been assumed with 
r = Q m h — 0.5. The set of overdensities corresponding to 
different R for a given position in the linear density field (this 
can also be thought of as corresponding to a test-particle or 
small mass element) is called a trajectory. 

Examples of trajectories created for the filters defined 
by Equations hi, ti & Ware presented in Fig. hi. These trajecto- 
ries were calcula ted us ing Langevin integration as described 
by Bond et al. (1991). First, a Monte-Carlo realisation of 
a Gaussian random field was created, integrated over con- 
tiguous spherical shells in fc-space. The amplitudes of the 
field integrated over non-overlapping shells are independent 
and, for thin shells from In k — > In k + 8 In k, are given by a 
Gaussian random variable with variance 



2 
°"shcll 



-^rk 3 P(k)6lnk. 



(5) 



For the trajectories presented in Sections W & N, the field was 
determined around a given position for 2 12 shells equally 
spaced in —7 < lnfc < 6. Trajectories were then created 
by summing over these shells, weighted by the filter win- 
dow function for 2 8 radii equally spaced in log a\ 4 . As a 
check that the density field is adequately sampled, results 
are presented in Sections kj & g for the sharp fc-space filter, 
for which the result is known analytically. Tests have also 
shown that increasing or decreasing the sampling does not 
significantly alter the results. A detailed analysis of the sam- 
pling required to calculate the backwards conditional mass 
function is given in Section M. 

In order to provide a complete description of the hi- 
erarchical build-up of mass, the epoch of collapse must be 
linked with the overdensity. In fact, the link between the 
critical overdensity S c and the collapse epoch is predicted 
to be independent of mass for the simple spherical top-hat 
collapse model (Percival, Miller & Peacock 2000) , the model 
adopted in standard PS theory. The mass of halo surround- 
ing a given position at time t is then given by the first, or 
largest mass, upcrossing of a horizontal barrier 8 — S c (t) by 
the trajectory. 



3 PARAMETERS IN PS THEORY 

This Section discusses the free parameters in Press- 
Schechter theory, and outlines the conventions used in this 
paper. 

The choice of filter was presented as a single free param- 
eter in Section El However, because PS theory only predicts 
an analytic form for the unconditional and conditional mass 
functions for the sharp fc-space filter, many authors have in 
the past used these analytic functions and have calculated 
ajyi using a different filter. The choice of filter used to calcu- 
late erf/ becomes an additional free parameter in PS theory 
(with the choice of filter to determine the mass function 
being a different free parameter). For instance, both Kauff- 
mann & White (1993) and Somerville & Kolatt (1999) use 
the spatial top-hat filter to calculate a\j, and the sharp fc- 
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Figure 1. Trajectories are plotted calculated by Langevin integration for three independent positions in a Gaussian random field with 
CDM power spectrum. For each position, an estimate of the density field was made integrated over 2 16 contiguous shells equally spaced 
in lnfc. Each Monte-Carlo realisation of the field was then filtered with sharp fc-space (solid line), top-hat (dashed line) and Gaussian 
(dot-dash line) filters at 2 12 radii, equally spaced in logcr^. The top-hat and Gaussian filters provide smoother trajectories than the 
sharp fc-space filter for the same field position although, as expected, the overall behaviour is similar for all filters. Trajectories are 
matched as described in Section H. The mass associated with the sharp fc-space trajectory shown in the centre panel will undergo a jump 
from a value corresponding to Ra to a value corresponding to Rg at a time when S c = Sj^g, the value along the dotted line. 
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Figure 2. The transitional weighting applied to the fc-shells for 
one step in a trajectory for the sharp fc-space (solid line), top-hat 
(dashed line) and Gaussian (dot-dash line) filters for a change in 
mass of between 10 and 5 X 10 Mq. The filters have been 
matched so that the transitional weighting peaks at the same 
value of lnfc, as described in Section H. 



space filter to determine the mass functions. Similarly, Lacey 
& Cole (1994), allowed the choice of filter used to calculate 
a\f to be different to the sharp fc-space filter used to deter- 
mine the mass function. The effect of this switch is equiv- 
alent to using the sharp fc-space filter throughout with a 
slightly modified power spectrum. In this paper, the same 
filter is used to determine both the mass function and <j\ { . 

Another parameter that is often allowed to vary is the 
critical overdensity for collapse, 8 C , with authors tuning this 
parameter to fit simulations (Lacey & Cole 1994), or using 
the ellipsoidal collapse model (Sheth, Mo & Tormen 2001). 
For simplicity, only the spherical top-hat collapse model is 
considered in this paper. 

In order to compare results calculated using different 



filters, the radii of the different filters need to be matched. 
In each step of a trajectory, the power spectrum is convolved 
with a transitional weighting function. Fig. ti shows an exam- 
ple of the transitional weighting applied by the three filters 
considered in this paper. Matching the radii of the differ- 
ent filters effectively matches the weightings applied to the 
power spectrum. A common method is to integrate each fil- 
ter over all real space to calculate its 'volume of influence' 
and the average mass enclosed. Equating volumes leads to 



Rks — 



(t) 



-1/3 



Rth 



0.4147? 



th: 



Rg — 



4tt\ 1 /3 



(f) 



(2*) 



1/2 



R th ~ 0.6437?^ 



(6) 



(7) 



For the differentiable filters it is also possible to match filter 
radii to second order in the expansion around fc = 0. This 
requires R g = l/y/ER th ^ 0.4477? th . 

Each point in a trajectory is calculated by integrating 
over fc-shells, and matching different filters simply alters the 
region of the power spectrum probed. It is possible to di- 
rectly match the peak value of fc probed at a given step 
in a trajectory. For the non-differentiable sharp fc-space fil- 
ter, fc m ax = 1/Rks- For the Gaussian filter the first solu- 
tion to the differential equation d 2 W / dRdk = is given by 
fcmax = V2/R g and for the top-hat filter solving this equa- 
tion numerically gives fc ma x — 3.342/ Rth- Combining these 
equations gives 



Rk 



: 0.299.R t 



R a ^ 0.4237? t , 



(8) 



The differential weighting applied in a particular trajectory 
'step' by each filter for this choice of matching filters is shown 
in Fig. g. 

In order to highlight the effect of filter choice on halo 
growth, it seems prudent to ensure that the same region 
of the power spectrum is probed at each trajectory step. 
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Figure 3. The halo mass function is presented for each filter as a function of v = S c /um for six values of S c , calculated from 10 6 
Monte-Carlo trajectories. The six values of <5 C correspond to redshifts of 0.0, 0.1, 0.5, 1.0, 2.0 and 4.0 for an Einstein-de Sitter cosmology. 
The prediction of PS theory for a sharp fc-space filter is also plotted (dotted line), and is a close match to the corresponding numerical 
data. The slight deviation is caused by differing sampling of the field and trajectory. Because of the filter shapes, such a difference will 
affect results for the top-hat & Gaussian filters less. Results for differ ent filters were matched as described in Sect ion H . The dashed line 
shows the fit to N-body simulation results of Sheth & Tormen ( f 99£ ), and the solid line the fit of Percival et al. (GOOOl). 



Because ? of this, where appropriate, the filter radii have been 
matched using Equation |8| Having adopted this convention, 
the halo mass M was defined as 

4tt . 



M 



-pRthi 



0) 



where p is the mean density of the universe. Because the 
mass variance, a\.j is dependent on filter choice (Equation W), 
comparing results as a function of a\ t implicitly matches 
filter radii. In this case, the relationship between the radii of 
different filters is dependent on the power spectrum. Filter 
radii are only matched using Equation H when the halo mass 
is required. 



4 THE HALO MASS FUNCTION 

The mass function f(M\t) dM, gives the distribution of mass 
in isolated halos at a given epoch, and is related to the num- 
ber density of halos by pf(M\t) dM = Mn(M) dM. The 
mass function can be calculated analytically for a sharp 
fc-space filter because the trajectories are Brown ian ran- 
dom walks and therefore have special symmetries (Press & 



Schechter 1974 Peacock fe Heavens 1990| [Bond et al. 199l[ ). 

(10) 



Defining 

S c 



the PS mass function can be written as 



/(ln^jt) dln^ = \l —ve 



2 /2 



d\nv. 



(11) 



Note that when the mass function is written in this way, its 
form is independent of epoch. It is not possible to obtain 
corresp onding analytic fo rmulae for top-hat or Gaussian fil- 
tering ( |Bond et al. 199 1| ). 

By creating an ensemble of 10 6 trajectories for each of 
the filters given byEquations nl & using the method de- 
scribed in Section bl the mass function has been estimated at 
six epochs. These data are presented in Fig. Q as a function 
of v. 

Recent work analysing the mass function predicted by 
numerical simulations has suggested that, although the mass 
function has a universal form when written as a fu nction 
of is, this form is altered from that of Equation hll (pheth 
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Figure 4. The distribution of halo creation times are plotted for six halo masses as a function of v = S c /rJM- These data were calculated 
from If) 6 Monte-Carlo traj ectories created as described in Section bl The lines are as in Fig. H, converted into c reatio n rates by multiplying 
by v (Percival, et al. 2000|). The Monte-Carlo data have been normalized to the model of Sheth & Tormen (1999). 



& Tormen 1999; |Percival, 
Instead, Sheth & Tormen (^999J) 
the N-body data was given by 



/<l'"'l'>-'l.M<=.L A /~(l +i ^ 



et al. 200C; Jenkins et al. 2001) 



found that a good fit to 



-i/ 2 /2 



din;/, 



(12) 



where v' — a x ' 2 u and a & p are parameters. A is determined 
by requiring that the probability density function is correctly 
normalized. Sheth & Tormen found best fit parameters a — 
0.707 and p = 0.3 for their simulations and group finding 
algorithm. The standard PS multiplicity fu nction has a — 
1, p = and A = 1/2. In Percival et al. (2000), a small 
tweak was required for the parameters in order to fit the 
data giving a = 0.774 & p = 0.274. These fits to the mass 
function are presented in Fig. H. 

The mass function calculated from trajectories deter- 
mined for the sharp fc-space has close to its predicted an- 
alytic behaviour given by Equation O, showing that the 
Monte-Carlo procedure works for this filter. The smoother 
trajectories calculated for the top-hat and Gaussian filters 
produce mass functions more like that observed in numerical 
simulations, with fewer high mass halos than predicted by 
Equation hll This match is not perfect, and the smoother 
filters predict more low mass halos than Equation u2. 



5 THE HALO CREATION RATE 

The halo creation rate f(t\M) dt for halos of mass M is given 
by the distribution of first upcrossing times of a barrier at 
the filter radius corresponding to M. The distribution of 
halo creation times determined for Gaussian, top-hat and 
sharp fc-space filters are plotted in Fig. H for the subset of 
Monte-Carlo trajectories that predict the creation of a halo 
with one of six masses. Each trajectory does not necessarily 
have a first upcrossing at a particular mass because the mass 
can undergo finite 'jumps' (see Section M). Consequently, 
the total number of creation events recorded from a fixed 
number of trajectories depends on the filter chosen. As the 
sharp fc-space filter provides a 'rough' trajectory compared 
to top-hat and Gaussian, there are more mass jumps and 
less creation events for this filter at a given final halo mass. 
In fact, it is only because of the discrete nature of these 
trajectories that any creation events at a particular mass 
are found. 

It follo ws fro m the Bayesian argument presented in Per- 
cival et al. (E000|) that 



i>/(lnf[t) dim; oc f(]nv\M) din v. 



(13) 



In the first of these probability density functions v is con- 
sidered to be a function of mass, while in the second it is 



© 0000 RAS, MNRAS 000, 000-000 



6 Will J. Percival 



] — i — i — i — i — | — i — i — i — i — p - 1 — i — i — i — | — i — r - 

Qi 



1 1 1 1 1 1 1 1 1 1 1 1 iDajj — i 1 1 1 r- 

o sharp A;— space 
A top— hat 
□ Gaussian 



l 1 1 1 1 1 1 1 1 1 1 1 1 r^n 1 1 r 



o sharp k— space 
A top-hat 
□ Gaussian 



ocsl 



o sharp fc-space 
& top-hat 

□ Gaussian A j<r 



CXfflD. 



S&, 



03 



,o° 



ivfect 



o o 

O O ' 



' v > 

\ ^ ■ 

4» O O 

3 °o. 'e 

7 a • 

::P 
o m 

.0 □ 



I" 



op 

I .... I 



A6 C =0.0169 

i i i I i i i i_ 



p"ni —■ - o 



caJ 



A6 =0.169 



■ ■ ■ ■ 



1111 



*~ 



11 




A6 =0.845 



-3 



-2 



-1 



-2 -1 

loQio^i^) 



Figure 5. The backwards conditional mass function calculated for PS theory with three filters is presented as a function of fi_>2 > for &\ 
equivalent to a mass of 10 14 Mq. All trajectories that predict a halo of mass 10 14 Mq at 82 > 1.69 were included in this analysis. The 
dotted line shows th e ana lytic prediction of PS theory with a sharp fc-space filter, and the dashed line shows the fit to the mass function 
of Sheth & Tormen (199E), presented as if its form is invariant to the change v —> fi— >2. 



considered to be a function of time. The model mass func- 
tions plotted in Fig. H have been converted to creation rates 
by multiplying by v and renormalising, and are compared 
with the creation rate calculated from the Monte-Carlo tra- 
jectories in Fig. 0. 

The results are similar to those of Fig. H, with the 
Monte-Carlo creation rate calculated for the sharp fc-space 
filter closely following the expected form. The creation rates 
calculated using the fit to the N-body simulations (Equa- 
tion R2) are closer to PS theory with the top-hat and Gaus- 
sian filters although, as for the mass function, the match is 
not perfect. Comparing Figs. H&W shows that the Bayesian 
link between the halo creation rate and mass function is 
closely followed by the Monte-Carlo data. 



6 THE BACKWARDS CONDITIONAL MASS 
FUNCTION 

6.1 PS theory 

The subset of trajectories that pass through a selected posi- 
tion determine the probability density functions for the past 
and future masses of that halo. As for the mass function, 
analytic formulae are possibl e for the sharp fc-sp ace filter 



(Bower 1991 



Bond et al. 1991 



Lacey fc Cole 1993), whereas 



for top-hat and Gaussian filters numerical techniques must 
be used. The conditional analogue of v is defined as 



Vi- 



81 — 82 



yfo 



(14) 



where a\ < a\ (M2 > Mi) and 82 < 8\ (£2 > ii), and the 
subscripts M and c have been omitted for simplicity. For the 
backwards conditional mass function, o\ and 82 (and 5i) 
known and ^1^2 is considered to be a function of a\. 
the sharp fc-space filter 

/(In 1/1^2 1 <5i , & , ct 2 ) d In ^1^2 



arc 
For 



^i^2e 



2/2 dlnj/i_ 2 



(15) 



which has the same form as Equation hll This is expected 
because the form of the trajectories is invariant to a linear 
shift in 8 or a 2 . 

For the top-hat and Gaussian filters, the equivalent dis- 
tribution can be calculated by analysing numerical realisa- 
tions of trajectories. In order that all upcrossings at small 
AS C = {82 — S\) were realized, the trajectories used for this 
analysis were sampled at 2 10 steps in R, equally spaced in 
—6 < log 10 (<T 2 — a 2 ) < 2. For the sharp fc-space filter, only 
a single fc-shell covering each interval in R is required. The 
other two filters have a broader shape and each step in a 
trajectory will depend on a larger region of fc-space. Two 
fc-shells within each interval in R and 2 7 shells either side 
of the minimum and maximum radii, equally spaced in lnfc 
were found to sufficiently sample fc-space. 

The backwards conditional mass functions resulting 
from these trajectories are presented in Fig. H for three val- 
ues of A8 C . In order to match the N-body results (for an 
Einstein-de Sitter cosmology) , only trajectories that predict 
a halo of mass 10 14 Mq at 81 > 1.69 have been counted. The 
smoother trajectories calculated using top-hat and Gaussian 
filters predict low mass progenitors, particularly for small 
AS C . 



6.2 Numerical Simulations 

Three numerical simulations, completed using the Hy- 



dra N-body, hydrodynamics code (Couchman, Thomas & 
Pearce 1995), have been analysed in order to estimate the 



backwards conditional mass function for halos of mass ~ 
10 14 Mq . Details of the cosmologies adopted for these sim- 
ulations are summarized in Table gj. Each simulation con- 
tained 128 3 dark matter particles, and particle positions 
were output at a large number of epochs (362 for FCDM, 
345 for OCDM & 499 for ACDM), separated by approx- 
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Figure 6. The backward conditional mass function calculated from three N-body simulations (See text for details). The left column 
shows the backwards conditional mass function from consecutive output from the simulation, while the centre column shows the backward 
conditional mass function for A<5 C ~ 0.169, and the right column A5 C ~ 0.845. The rows show results from different simulations. Lines 
are plotted as in Fig. H. 



model 



TCDM 


1 





0.25 


0.64 


0.5 


OCDM 


0.3 





0.15 


0.85 


0.5 


ACDM 


0.3 


0.7 


0.15 


0.85 


0.5 



Table 1. Table showing the parameters of the cosmological mod- 
els adopted in the three N-body simulations. 



imately equal intervals in time. The box size chosen was 
100/i -1 Mpc for all three simulations that gave a particle 



mass of 2.6 x 10 11 M for TCDM and 7.9 x 10 10 M for the 
other two simulations. Groups of particles were found using 
a standard friends-of-friends algorithm with linking length 
set to b — 0.2 times the mean interparticle separation. 

Fig. H shows the conditional mass function calculated 
for different values of A8 C . Halos of between 200 and 1500 
particles were selected from every output time, correspond- 
ing to masses 1.6 x 10 13 M < M < 3.9 x 10 14 M , and the 
particles were traced back to determine the mass of the halo 
they resided in at the output time nearest the required A5 C . 
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Only halos containing at least 40 particles, and with over 
50 % of their particles in common with the final halo were 
considered to be progenitors. The choice of 50 % is reason- 
ably arbitrary, and does not affect the progenitor distribu- 
tion significantly. a\[ was calculated for these halos using 
the top-hat filter. We have tried using Gaussian and sharp 
fc-space filters to calculate a 2 M and find that the result is not 
significantly affected by this choice. 

The values of A<5 C adopted correspond to Az ~ 0.01, 
Az = 0.1 and Az = 0.5 for an Einstein-de Sitter cosmology. 
The smallest of these intervals corresponds to consecutive 
output from the simulations. 

Comparing Fig. and Fig. H provides weak evidence 
that the build-up of halos in PS theory is most like that 
predicted by the sharp fc-space filter. This is discussed fur- 
ther in Section [j. 



7 THE COSMOLOGICAL DISTRIBUTION OF 
MERGERS 

In order to link halo growth with the observed evolution of 
the baryonic component of galaxies, it is important to dis- 
tinguish between major mergers, which can disrupt struc- 
ture creating observable signatures, and the slow accretion 
of matter. Where trajectories dip and rise again (as the fil- 
ter radius decreases), the mass associated with the trajec- 
tory undergoes a jump. Such an event is shown in Fig. hi 
for a trajectory calculated with a sharp fc-space filter. These 
jumps are often associated with mergers, and it is expected 
that close congregations of excess overdensity in the linear 
density field that initially collapse as single entities and sub- 
sequently merge together, will create such a dip and rise as 
the filter increases in size and the number of enclosed clumps 
increases. For PS theory with a sharp fc-space filter, jumps 
are the only mechanism for gaining mass as the trajectories 
are not differentiable. 

The halo creation events analysed in this Section can be 
thought of as being instantaneous. However, numerical lim- 
itations are placed on how well we can follow these events 
because the trajectories are discretized in mass and the sim- 
ulations are discretized in both mass and time. For the sim- 
ulations, simultaneous outputs from the simulations were 
used to determine whether a halo had just been created, and 
whether the creation event was a merger. For the trajecto- 
ries, the behaviour of a trajectory before it passed through 
a point was analysed to determine whether it was a first up- 
crossing (a creation event), and the subsequent behaviour 
was analysed to determine the form of creation. Both the 
trajectories and the simulations were sufficiently sampled 
that this numerical problem should not affect the conclu- 
sions drawn from this work. Defining the merger fraction, 
fm*rg(Mf,Zf), by 



f mel g{M f ,z f ) = 



No. mergers 



No. new halos 



(16) 



where Zf is the final redshift of the halo and Mf is the final 
mass of the halo, we can write the normalized proportion of 
new halos created by merger as 



-$ — &— 



o sharp A:— space 
a top— hat 
a Gaussian 

i I 



i i 



f*'(Mf,z f ) = 



fmci g (Mf,Zf) 
/mergCM/) 



(17) 



2 4 6 8 

redshift 

Figure 7. The normalized proportion of halo creation events that 
are a result of merger, measured from a set of 10 6 trajectories, 
calculated for each of the three filters described by Equations M - 
M. Here, a merger was defined to be a mass jump of at least a 
quarter of the final mass, with the trajectory after the jump also 
corresponding to a halo with mass over a quarter of the final 
value. 



where / mer g(Af/) is the average of f ulels (Mf, Zf) over all red- 
shift. These data were combined over all final halo masses of 
interest, weighted by the standard Binomial variance treat- 
ing the number of creation events found for each Mf and 
Zf as the number of 'experiments' performed. The result- 
ing normalized merger fraction, a function of redshift only, 
is denoted by fi. This procedure is designed to allow the 
merger fraction to be a function of mass and to highlight 
any cosmological changes in the overall merger rate. 

Fig. u\ shows the cosmological variation in fj,, calculated 
from PS trajectories, for mergers defined such that both 
the mass jump and the progenitor correspond to at least a 
quarter of the mass of the final halo. Results presented in 
this Section are not sensitive to this choice. The numbers of 
mergers and halo creation events were determined from a set 
of 10 6 trajectories, calculated as described in Section H for 
each of the three filters given by Equations hi- pi The discrete 
nature of the trajectories provides the required lower limit of 
the size of possible accretion event and stops the number of 
halo creation events in any time interval becoming infinite. 
Data have been combined for all creation events that result 
in a halo of mass 4 x 10 10 M© < M < 10 13 M , in such 
a way that the results should be directly comparable with 
those calculated from N-body simulations. 

The Markovian nature of the trajectories calculated us- 
ing a sharp fc-space filter means that the probability that a 
halo was created by merger is independent of creation time: 
the normalized merger probability is flat. However, the tra- 
jectories calculated for the other filters are non-Markovian 
and the probability of merger is dependent on cosmolog- 
ical epoch. There is a general trend for these filters that 
when a halo is created at early times it's more likely to have 
been created via the slow build-up of matter rather than by 
merger. 
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Figure 8. The normalized proportion of halo creation events that 
are a result of merger, measured from three N-body simulations 
from which groups of particles were determined using a standard 
friends-of-friends group finder with b = 0.2. Halo creation was 
defined to have occurred for a given halo at a particular time if 
at least half of the particles had not previously been contained 
within a halo with equal or greater mass. Only halos comprising 
at least 45 particles were included in this analysis. A creation 
event was said to have been a merger if there were two possible 
progenitors at the previous output time from the simulation, each 
with mass greater than a quarter of the final halo mass. The data 
are not strongly dependent on this choice, and are consistent with 
no evolution in the proportion of merger events as a means of halo 
creation. 



It is interesting to compare this behaviour with that 
observed in N-body simulations. Here, the discrete nature 
of the simulations stops the number of halo creation events 
within a time interval from tending to infinity. Halo creation 
was defined to have occurred for a given halo at a given 
time if at least half of the particles had not previously been 
contained within a halo with equal or greater mass. Only 
halos comprising at least 45 particles were included in this 
analysis. A creation event was said to have been a merger if 
there were two possible progenitors at the previous output, 
each with mass greater than a quarter of the final halo mass. 
Using the analysis described above for combining data from 
halos of different mass, these data are presented in Fig. H. 

Fig. ra suggests that the distribution of mergers follows 
the PS model with a sharp fc-space filter, with no evidence 
for a cosmological trend as predicted by PS theory with 
either the top-hat or Gaussian filters. PS theory with the 
sharp fc-space filter would therefore seem to provide the best 
model for the build-up of clumps of dark matter determined 
using the FOF algorithm. 



8 THE TIE-IN WITH OBSERVED GALAXY 
EVOLUTION 

A variety of different techniques have been used to deter- 
mine the amount of evolution in the field galaxy merger 
rate as a function of redshift: recent studies have used ob- 



galaxy morphologies ( Lavery et al. 1996| ). See Abraham et al. 
(19981) for a review. There is now strong evidence that the 
number of interacting systems grows rapidly with look-back 
time. The cosmological change in the proportion of galax- 
ies involved in mergers is usually fitted by a power-law of 
the form (1 + z) n , and results are parametrized by the best 
fit value of n. A recent study using follow-up HST imaging 
of a sample of 185 galaxies with measured redshifts from 
the CFRS and Autofib-LDSS surveys has provided the mos t 
convincing measurement of n to date (Le Fevre et al. 2000). 
This w ork extends the pair-count analysis of Patton et al. 
(1997) from z — 0.33 to Z — 0.9 and reveals a consistent 
amount of evolution between such studies and results from 
the morphological classification of galaxies. A best fit value 
of n — 2.7 ± 0.6 is determined after correction for biases 
attributed to selection effects: the combination of ground- 
based and HST observations means that the sample is bi- 
ased towards finding close pairs that are thought to be single 
galaxies in the magnitude limited ground-based sample, but 
which can be resolved into separate galaxies by HST. 

These observational studies select galaxies based on 
properties of the baryonic galaxy component. In order to 
model such observations, the predicted evolution in the halo 
merger rate has to be combined with the time-scale over 
which the tracers of halo merger are visible, including any 
time difference between the merger of the baryonic compo- 
nents and the halos caused by the dynamical friction time- 
scale. The dynamical friction time-scale describes how long 
it takes for a galaxy (or progenitor) to spiral to the centre of 
a halo. If this takes a long time compared to the cosmological 
time-scale of interest, it would affect the relative distribution 
of baryonic and halo mergers. However, by concentrating on 
galaxies in small groups or field galaxies at relatively low 
redshif ts, the dynamical friction t imescale is expected to be 
small ( Barnes fc Hernquist 1992 ). For such studies, or by 
using galaxy counts rather than morphologies, the observed 
evolution is expected to be driven by the proportion of halos 
that have undergone recent merger. 

Percival et al. (200C) showed that multiplying the mass 
function by dS c /dt leads to the joint distribution of halo cre- 
ation events in mass and time. In Section tn the distribution 
of mergers in numerical simulations was shown to follow that 
predicted by PS theory with a sharp fc-space filter, and to 
have the same form as the distribution of creation events. 
The proportion of halos that have just undergone merger is 
given by the creation rate divided by the number of halos, 
given by the standard PS mass function (Equation |ll| ). The 
resulting cosmological change is proportional to dS c /dt, and 
only depends on the Friedmann equation. A meth od for cal- 
culating dS c /dt was given by Percival et al. (2000). Because 
PS theory does not include sub-halos that form part of a 
larger collapsed system, this is only valid for field galaxies. 
This will also be the evolution of any combination of halos of 
different mass, provided the combination is independent of 
redshift. For the three cosmologies parametrized in Table fil 
approximating dS c /dt by (1 + z) n gives n = 2.5 for SCDM, 
n = 2.1 for OCDM and n = 2.3 for ACDM. These values 
arc in agreement with the best fit value of Le Fevre et al. 
(P000|) , n — 2.7 ± 0.6 within the errors and uncertainties. 
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Further observations to provide more accurate observational 
data could, in principle, provide an interesting mechanism 
for constraining the best-fit cosmology. 



9 CONCLUSIONS 

The build-up of halos predicted by PS theory has been 
shown to be strongly dependent on the filter used to smooth 
the density field. None of the filters provides a perfect 
match for the build-up of structure in the numerical simula- 
t ions presented in this paper. The simulatio n mass function 
(Sheth & Tormen 1999; Percival, el al. 2000) is most similar 
with PS theory with a top-hat filter, as is the halo creation 
rate. 

The top-hat and Gaussian filters predict mass growth 
due to continuous accretion with the occasional merger, 
while the sharp fc space filter predicts mass growth com- 
pletely due to mergers. Because of this, the predicted back- 
wards conditional mass function is strongly dependent on 
the choice of filter. Analysing numerical simulations sug- 
gests a backwards conditional mass function most like that 
of PS theory with the sharp fc-space filter, although again, 
the match is not perfect. The resolution of the simula- 
tions means that it is difficult to resolve mergers with small 
enough halos. However, there is weak evidence that the sim- 
ulations contain more mergers with small mass halos than 
are predicted by the smoother trajectories associated with 
the top-hat and Gaussian filters. 

For the sharp fc-space filter, the cosmological distribu- 
tion of halo mergers is independent of merger parameters. 
This result follows from the Markovian nature of the tra- 
jectories produced. For the other filters, which predict non- 
Markovian trajectories, the probability of merger is shown 
to be dependent on the halo creation epoch. No evolution 
is seen in the probability that a given halo was created by 
merger measured from N-body simulations, and PS theory 
with a sharp fc-space therefore provides the best fit to these 
data. This provides the most convincing evidence that PS 
theory with a sharp fc-space filter provides the best model of 
halo growth. In this work, a standard FOF group finder was 
used. An interesting test would be to analyse the observed 
build-up of halos using a number of different algorithms to 
find halos. However, this is not expected to significantly alter 
the conclusions because the build-up of mass in the simula- 
tions is independent of the group finder. 

When combined with the link between the mass func- 
tion and halo creation rate described in Percival el al. 



( 2000), this leads to the interesting result that the cosmolog- 
ical change in the proportion of halos that have undergone 
recent merger is predicted to be independent of power spec- 
trum, halo mass, mass function, and to only depend on the 
Friedmann cosmology through d5 c /dt. In Section H it was ar- 
gued that the field galaxy merger rate should be driven by 
the halo merger rate, and that observations of field galax- 
ies estimate the cosmological evolution in the halo merger 
rate. In Section kI recent observational data were shown to 
be consistent with the presently favoured ACDM cosmology, 
although they are insufficiently accurate to usefully distin- 
guish between cosmological parameters. It is hoped that fu- 
ture studies of this evolution will provide sufficiently accu- 
rate results that this picture will change. 
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